Cold Attractive Spin Polarized Fermi Lattice Gases and the Doped Positive U 

Hubbard Model 
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Experiments on polarized fermion gases performed by trapping ultracold atoms in optical lattices, 
allow the study of an attractive Hubbard model for which the strength of the on site interaction 
is tuned by means of a Feshbach resonance. Using a well-known particle-hole transformation we 
discuss how results obtained for this system can be reinterpreted in the context of a doped repulsive 
Hubbard model. In particular we show that the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
corresponds to the striped state of the two-dimensional doped positive U Hubbard model. We then 
use the results of numerical studies of the striped state to relate the periodicity of the FFLO state 
to the spin polarization. We also comment on the relationship of the d x 2_ y 2 superconducting phase 
of the doped 2D repulsive Hubbard model to a d-wave spin density wave state for the attractive 
case. 

PACS numbers: 71.10.Fd, 03.75.Ss, 74.25.Ha 



Using standing wave laser light fields and optical Fes- 
hbach resonances, ultracold atomic gases can be used to 
realize a variety of Hubbard like models*^ Technically if 
one wants a strong interaction, it appears easier to cre- 
ate an ultracold fermi gas (e.g., 6 Li) with an attractive 
interaction and experiments on polarized fermion gases 
with attractive interactions have been carried out to look 
for exotic superfluid states.— 

One state of particular interest is the FFLO state 
which arises from pairing across the spin-split Fermi 
surface of a spin polarized systenijL^ In this state the 
Cooper pairs have a finite center of mass momentum 
leading to spatial oscillations of the pair field order 
parameter and the spin polarization. Here, making 
use of density-matrix-renormalization-group (DMRG) 
results^ - for the doped, positive U Hubbard model and a 
well-known particle-hole transformation ; 11 : 12 ' 13 : 14 we dis- 
cuss the properties of the FFLO state for the case of a 
half-filled spin polarized negative U Hubbard model on 
a two-dimensional lattice. We also comment on the re- 
lationship of this to the question of whether the doped 
positive U Hubbard model may have a d x i_ y i supercon- 
ducting ground state. 

The Hamiltonian for a half-filled negative U Hubbard 
model in an external Zeeman field, H, can be written as 



« = -*£(4 



nil)' 



(1) 



Here t is a nearest neighbor one-electron hopping, U is 
positive so the onsite interaction is — f7nji-nji and h = 
gHoH/2. The operator c\ s creates a fermion with spin s 
on site i, the sum (ij) is over nearest neighbor sites of 



a square lattice and n, s — c\ s Ci S is the number operator. 
In the case of the ultracold gases, a spin polarization 
p = (JVf — N±)/N is achieved by loading more atoms in 
the "up" pseudo-spin hyperfine state than in the "down" 
state. 

In the mean-field FFLO state, the Cooper pairs have 
a finite center of mass momentum which, for a two di- 
mensional tight binding band structure, lays along the 
(1,0) or the (1,1) direction depending upon the ratio of 
\U\/t. For the case in which it lays along the (1,0) di- 
rection, the order parameter exhibits a one-dimensional 
oscillation along the x direction 



A(4) = Re(c^c ei ) = -A cos(g x 4). 



(2) 



Here we have chosen the phase so that A(0) is negative. 
The spin polarization also varies in space with 



n s {t x ) = (n £T —nti) = p-m cos(2q x £ x ). (3) 

Here cos(q x £ x ) and cos(2q x £ x ) give the leading weak cou- 
pling harmonic behavior. The spatial variation of A (£ x ) 
and n 3 (£ x ) are schematically illustrated in Fig. [TJ At 
stronger coupling higher harmonics can enter giving a 
more localized behavior. 

Since the possibility of an FFLO state was originally 
proposed,^ there has been great interest in determining 
whether it exists and exploring its properties. In par- 
ticular, with the possibility of observing such a state in 
ultracold atomic gases, there have been a number of the- 
oretical calculations . 15 : 16 i 17 : 18 : 19 : 20 : 21 These calculations 
have been based upon a mean field description. Here we 
are interested in the case of a two-dimensional half-filled 
lattice gas where it is important to treat strong coupling 
effects. One would like to know if the FFLO state sur- 
vives beyond the mean field approximation, whether the 
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stripes run vertically (or horizontally) as in Eqs. ([2]) and 
© or diagonally, and whether the periodicity of the spa- 
tial variation which determines q x remains at its mean 
field value. 

The Hamiltonian for the negative U Hubbard model 
in a Zeeman field can be transformed under a unitary 
transformation to a doped positive U Hubbard model. 
This transformation, introduced by Emery, 11 has been a 
staple in quantum Monte Carlo work where it was used to 
show that the fermion sign problem is absent for the half- 
filled Hubbard model with a near- neighbor hopping. 12 
It has also been used, as we will review, to provide a 
map between the half-filled negative U Hubbard model 
in a Zeeman field and the doped positive U Hubbard 
model . 12 ' 13 ! 14 ' 22 Under this transformation we will see 
that the FFLO state becomes the familiar striped, charge 
density, and 7r-phase shifted antiferromagnetic state of 
the doped positive U Hubbard model. 

For a unitary transformation in which 



the Hamiltonian given in Eq. |1} becomes 



(4) 



n = -tJ2 «W + 4-*-) + 2 ^ ( " lT + Utl ~ 1)5 



Mean field studies 1 ^ 3 -^^^! 2 ! of the doped, two- 
dimensional positive U Hubbard model with a nearest- 
neighbor hopping find that the domain-walls run along 
the x or y-directions for \U\/t < 3.6. For values of 
3.6 < \U\/t < 8 the domain-walls run diagonally and 
for \U\/t greater than 8 the domain- walls are found to 
be unstable with respect to the formation of magnetic 
polarons^ The mean field domain walls are found to 
contain one hole per site along the wall and one says 
that the domain wall is filled with holes. This means 
that the spacing of the vertical (or horizontal) charged 
domain walls is equal to giving the mean field result 



q x = K-p. 



(9) 



(5) 



for such walls. A schematic illustration of the mean field 
results for p = 0.125 is shown in Fig. QJi. 
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with [a = —h. That is, the half-filled negative U Hubbard 
model with spin polarization p is mapped into a positive 
U Hubbard model with a site filling (n) = 1 — p. Under 
this transformation Eq. ^ becomes 



(-iy* +t »(M x (£ x J y )) = -A cos fe 4 



(6) 



with M x (£) = (d\^dii + d\yd^)/2. The system has ro- 
tational symmetry so that Eq. © implies that the stag- 
gered spin order is modulated by cosq x £ x . 



m st ag(£x) = ^ («-l(4,^y) ~ ™l(4,^,)} 

= — A cos q x t x . 



(7) 



Similarly, the spin polarization, Eq. ([3]), transforms to 
the hole density giving 

h(£ x ) = 1 - (n(£ x ,£ y )) =p-m cos(2q x £ x ). (8) 

Thus the FFLO state maps to the striped state of 
the doped positive U Hubbard model in which charged 
stripes separate 7r-phase shifted antiferromagnetic re- 
gions. 



FIG. 1: (Color on line) Schematic results for the spatial varia- 
tion of the s-wave order parameter A(£ x ) (solid) and the spin 
polarization n s (£ x ) (dashed) in an FFLO state with p — 0.125. 
(a) The mean field result with a period for A(£ x ) of 2/p = 16 
and (b) with a period 1/p = 8 expected from DMRG calcu- 
lations on long stripes. 



Now these are mean field results and at larger values of 
U one is dealing with a strongly correlated system. Thus 
one might wonder whether the FFLO state remains the 
ground state. This is a question that has also played 
a central role in the high T c cuprate problem. There 
one would like to know if the ground state of the doped 
Hubbard model is striped or possibly a d x 2_ y 2 super- 
conductor, and what is the nature of the interplay be- 
tween stripes and pairing. 29 Present DMRG calculations 
on 6-leg Hubbard ladders find a striped ground state^i 1 ^, 
while dynamic cluster Monte Carlo calculations 3 ^ on pe- 
riodic lattices find evidence of a d x 2_ y 2 superconducting 
state. Whether the difference of lattice aspect ratios and 
boundary conditions or subtle numerical biases are re- 
sponsible for this disagreement is not known. It does, 
however, appear that these two states are very close in en- 
ergy. Here we will use the DMRG results for the striped 
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state to discuss the nature of the FFLO state and then 
conclude by noting what would happen for the d x i_ y i 
state. 

The DMRG calculations have been carried out on 6- 
leg doped positive U Hubbard ladders. Periodic bound- 
ary conditions have been imposed in the 6-leg direction 
and open boundary conditions are used in the long in- 
direction. In Ref. 10, the results have been extrapo- 
lated to infinite length ladders. These numerical results 
find that the ground state exhibits the striped structure 
shown in Fig. [2] for a 6 x 21 ladder with U/t = 12 and 
12 holes corresponding to a doping p = 2/21. The figure 
shows the hole density 

1 6 

M4) = r X)(l-<n(4,*v)», (10) 
e y =i 

and the staggered spin density 
i ! 

m stag (4) = gj] (-r/*+^(n T (4A) -ni(t x J y )}, 

ly=l 

(11) 

versus i x . For this doping, extrapolated DMRG results^ 
show that the ground state is striped for U/t > 3. For 
a 6-leg tube, the stripes are found to contain 4 holes 
(2/3 filled stripes) so that the charge stripe spacing is 
2/(3j?) = 7 with (n) = 1 -p. 




FIG. 2: (Color online). This figure shows DMRG results 
(G. Hager et al.— ) for the hole h(£ x ) (dashed red) and the 
staggered spin m B ta a {(. x ) (solid blue) densities along the leg 
direction for a 21 x 6 Hubbard ladder with 12 holes and U/t = 
12. As discussed in the text, h(£ x ) corresponds to the spin 
polarization n 3 (£ x ) and m a t ag (£x) corresponds to the s-wave 
pairfield order parameter A(£ x ) of the FFLO state. 

Calculations on an 8-lcg t — J ladder find that the 
stripes are half-filled and there is evidence that this is 
the preferred filling. 3 ! These calculations show that the 
holes that make up the stripe exhibit short range pairing 
correlations 3 ^ in which the sign of the singlet pair field 



associated with two sites connected by a hopping per- 
pendicular to the stripe is opposite to that of the case in 
which the hopping is parallel to the stripe. These d-wave 
like pairing correlations of the holes in the stripe are not 
taken into account in the mean field calculations which 
find filled stripes. For half-filled stripes, the spacing be- 
tween the charged stripes is l/2p so that the numerical 
calculations on the longer stripes imply that q x for the 
FFLO state will be twice that found in the mean field 
calculations. 

q x = 2vrp. (12) 

In this case, the FFLO state will have half the spa- 
tial period of the mean field solution, as illustration in 
Fig. [lb. 

Here we have made use of DMRG calculations which 
were carried out on systems with periodic boundary con- 
ditions in one direction and open boundary conditions 
in the other. In cold atoms experiments, one is gener- 
ally dealing with an underlying optical lattice in a slowly 
varying trap potential which acts as a spatially depen- 
dent chemical potential. This leads to soft boundaries 
and the simultaneous coexistence of spatially separated 
phases. The details of this depend upon the total filling 
as well as the curvature of the confining trap potential. 33 
In principle, the DMRG method is well suited to treating 
such soft boundaries but here we have focused our dis- 
cussion on the idealized case of a two-dimensional lattice 
which has one fermion per site with an attractive Fes- 
hbach tuned interaction and a spin polarization p. We 
have shown that the FFLO state of a 2D half-filled at- 
tractive spin polarized Hubbard model is a unitary trans- 
form of the striped state of the doped repulsive Hubbard 
model. Thus an experimental observation of the FFLO 
state in the above mentioned cold atom set-up would 
imply that the 2D doped repulsive Hubbard model is 
striped. Note that we have not shown that the striped 
state is the ground state of the two-dimensional doped 
repulsive U Hubbard model. This remains an open ques- 
tion. We have simply shown that if the striped state were 
the ground state then this would imply that the FFLO 
state would be found in the optical lattice experiments 
and that its periodicity would be given by Eq. (TT2"|) . We 
believe that the factor of two increase in q x represents 
an important change from the mean field result and re- 
flects the tendency towards d x 2_ y 2 pairing^ Recent ex- 
periments also suggest that local pairing correlations are 
present on the cuprate stripesi« 

Finding evidence for an FFLO state in a cold attractive 
spin polarized fermi lattice gas would provide evidence 
that the low temperature phase of the doped positive 
U Hubbard model is striped. If on the other hand, the 
d x 2_ y 2 superconducting state describes the low temper- 
ature phase of the repulsive Hubbard model, then the 
transformation, Eq. Q, tells us that the cold atom sys- 
tem will exhibit a "d-spin density wave" ground state 
characterized by an order parameter 
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^2(cosk x - cosk y ){c\ +m c ki ) (13) 
k 

with Q = (n,Tt). Thus experiments on half-filled, spin 
polarized cold fermi lattice gases with attractive on site 
interactions and the search for the FFLO state have im- 
portant implications for the properties of the doped re- 
pulsive U Hubbard model. Likewise, our present under- 
standing of the repulsive Hubbard model can provide in- 
formation relevant to experiments on the spin polarized 
attractive fermi gas. 
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